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Abstract 



Context. Surface abundances of OB stars provide important information on mixing processes in stellar interiors and on evolutionary 
status. Performing abundance determination from selected C, N and O lines in the optical spectra of OB stars is thus necessary to 
understand the evolution of massive stars. 

Aims. In this study, we investigate the formation of C m 4647-5 1-50 and C in 5696 in the atmosphere of O stars to see if they can be 
reliably used for abundance determinations. 

Methods. We use atmosphere models computed with the code CMFGEN. The key physical ingredients explaining the formation of 
the C m lines are extracted from comparisons of models with different stellar parameters and through examining rates controlling the 
level populations. 

Results. The strength of C m 5696 critically depends on UV C m lines at 386, 574 and 884 A. These lines control the C m 5696 upper 
and lower level population. C in 884 plays a key role in late O stars where it drains the lower level of C in 5696. C m 386 and C in 574 
are more important at early spectral types. The overlap of these UV lines with Feiv 386.262, Feiv 574.232 and S v 884.531 influences 
the radiative transfer at 386, 574 and 884 A, and consequently affects the strength of C m 5696. C m 4650 is mainly controlled by 
the Cm 538 line which acts as a drain on its lower level. Feiv 538.057 interacts with Cm 538 and has an impact on the Cm 4650 
profile. Low temperature dielectronic recombinations have a negligible effect on the line profiles. The main effect of the wind on the 
C m lines is an increase of the emission as mass loss rate increases, due to a shift of the formation region towards high velocity. Such 
a shift also implies desaturation of Cm 386, Cm 538, Cm 574 and Cm 884 which in turn modify the Cm 4650 and Cm 5696 line 
profiles. Metallicity effects are complex and result from an interplay between ionization, metallic line interaction with UV C m lines, 
and wind effects. 

Conclusions. Given our current understanding of the stellar and wind properties of O stars, and in view of the present results, the 
determination of accurate carbon abundances from Cm 4647-51-50 and Cm 5696 is an extremely challenging task. Uncertainties 
lower than a factor of two on C/H determinations based only on these two sets of lines should be regarder as highly doubtful. Our 
results provide a possible explanation of the variability of C in 4650 in Of?p stars. 

Key words. Stars: massive - Stars: atmospheres - Stars: fundamental parameters - Stars: abundances 
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1. Introduction 
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Massive stars turn hydrogen into helium through the CNO cycle 
in their core. In this process, nitrogen is produced at the expense 
of carbon and oxygen. Mixing processes subsequently transport 
the products of nucleosynthesis from the core to the outer lay- 
ers and surface of the star. The spectrum of a massive star thus 
contains the imprints of nuclear reactions taking place in the in- 
terior of the star. The older a star, the larger its nitrogen content 
and the lower its carbon abundance (provided the star is of type 
O and not yet an evolved object such as Wolf-Rayet (W-R), red 
supergiant or Luminous Blue Variable). The surface chemical 
composition depends on the star's rotation rate, metallicity, and 
mass. A fast rotation, for example, favors the transport of chem- 
ical species and thus quicker en richment of the photosph ere by 
material synthesized in the core dMaeder & MevnefeOOOh . Low 
metallicity stars burn hydrogen more efficiently, as does a 60 
M0 stars compared to a 20 M0 star. Chemical surface abun- 
dances are thus a key to understand the physical processes con- 
trolling the evolution of massive stars. It was surface abun dances 
that provided key clues that W-R stars were evolved (e.g JSmithl 



1973b ISmifh & Wilrislll982l). Recent results dHunter et al.l|2008h 



Przvbilla et al ■l l2010fclBrott et al Jl20 1 it iMartins et al.ll20 1 2^ have 
shown that strong constraints on the physics of massive stars 
could be obtained by spectroscopic analysis and abundance stud- 
ies. 

Nevertheless, the determination of surface abundances of 
O stars is challenging. Non-LTE atmosphere models including 
line-blanketing and winds are required, and there is a rarity of 
metallic lines in O star spectra. The UV range contains a few 
strong resonance lines (Cm 1176, Nv 1240, Civ 1550, Niv 
1720) but they are mainly formed in the wind part of the atmo- 
sphere where numerous effects (X-rays, clumping) render their 
modelling difficult. Photospheric lines are found in the optical 
range, but, here again, only a few lines are available. 

iRivero Gonzalez et al.l (1201 ll) have revisited the formation 
mechanism of N m 4634-4640 and have used several N n, N m 
and N I V lines to derive abun dances of O stars in the Magellanic 
Clouds. IMartins et alJ (1201 2l) studied the effect of magnetism on 
the surface abundance of O stars, and based their nitrogen abun- 
dance measurements on N 11, N m and N iv optical lines. Nitrogen 
is a crucial indicator of chemical enrichment. But more precise 
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constraints can be brought to stellar evolution if the ratio N/C 
can be derived because N/C is less sensitive to the initial metal 
content than N/H. For instance the ratio of N/C in the Galaxy to 
N/C in the SMC is less that 3, while N/H in the Galaxy is 13 
times larger than N/H in the S MC. Thus N/C, espec ially when 
combined to the N/O ratio (e.g. iPrzvbilla et al.ll20loh . provides 
stronger constraints on the enrichment and mixing history of 
CNO material at a given metallicity. 

Unfortunately, carbon lines are even less numerous than ni- 
trogen lines in O stars spectra. In the optical range, the strongest 
lines are Cm 4647-50-51 (hereafter C m 4650), Cm 5696, and 
Civ 5801-12. Cm 4650 is observed in absorption in most O 
stars. It shows emission only in the earliest O supergiants. In 
early O dwarf s, it is almost absent . A peculiar class of objects - 
the Ofc stars dWalborn et al .1120101) - displays always C m 4650 
emission. Stars of the Ofc class encompass Of?p stars which, 
among other characteristics, have a stron g but varying Cm 
4650 emission. Many of them are magnetic ( Donati et all 2006t 
iHubrig et aUl2008t iMartins et al.ll2010t IWade et al.ll201 lab . Ofc 



stars are found mainly at spe ctral type 05 but at all luminosity 
classes dWalborn et al.l l2010b . So far, no explanation has been 
given for the appearance of such C m 4650 emission. Its analysis 
is rendered difficult by the presence of the neighbouring strong 
Nm 4634-4640 lines. Civ 5801-12 is present in most O stars, 
but is often not observed. In late-type stars it is in absorption, 
while in early type stars it can exhib it a complex profile with 
broad emission wings and absorption dContill 19741 : iBouret et al.l 
120121) . 

C m 5696 is present in most O stars except the earliest ones 
(IWalbornlll980l) . It is in emission in giants and supergiants, and 
in most dwarfs. The e mission is o n average stronger in super- 
giants than in dwarfs dContill 19741) . Contrary to Cm 4650, Cm 
5696 is located in a spectral region free of other lines and hence 
its profile is well defined. Spectroscopic analysis have generally 
not considered this line for carbon abundance determination - 
atmosphere models often have difficulties reproducing the ob- 
served profiles, even if the H, He i, He n and UV CNO lines are 
well fitted. This indicates that the formation of that line is com- 
plex and pr obably depends on m issing or improperly accounted 
for physics. iNussbaumerl (Il97ll) found that the C m 5696 emis- 
sion was due to a depopulation of the lower level by C m 884. 
Hence, a correct description of the radiative transfer close to 884 
A is required. Since the models used by INussbaumerl (1 1 97 ll) did 
not include a proper line-blanketing and non-LTE treatment, we 
revisit the formation of C m 5696. We also investigate the forma- 
tion of C m 4650. 

In Section [2] we present the method and models we have 
used. Sect. [3] describes the formation mechanism of C m 5696 
and Sect. |4] focuses on C m 4650. The wind effects are gathered 
in Sect. [5] We summarize our results in Sect. [6] 

2. Model description 

Our investigation of the formation mechanisms of Cm 4650 
and Cm 5696 relied on the computation of model atmo- 
sphere and the analysis of the processes populating the lev- 
els involve d in the Cm transitio ns. We have used the code 
CMFGEN (iHillier & Millerll 19981) to compute atmosphere mod- 
els. CMFGEN produces non-LTE models including winds and 
line-blanketing. The radiative transfer is solved in the comoving 
frame and the geometry is spherical. The rate equation and the 
radiative transfer equation are solved simultaneously and itera- 
tively until convergence of the level populations. A super-level 
formalism is used to reduce the number of levels involved in the 



computations (to reduce the amount of memory required). Table 
lA.ll provides the number of levels used in our computations. Our 
models include H, He, C, N O, Ne, Si, S, Ai Ca F e, Ni. The ref- 
erence solar abundances of iGrevesse et al ] d2007l) were adopted. 
The radiative equilibrium condition sets the temperature struc- 
ture and the density structure in the inner atmosphere is obtained 
by iteration. Once the level populations are (roughly) converged, 
the radiative force is computed and the momentum equation is 
solved to yield the density structure; the radiative field and level 
populations are then re-computed until full convergence of the 
model is achieved. A couple of these hydrodynamical iterations 
are performed. In the out er atmosphere, the vel ocity law is given 
by a classical f3 law (e.g. lKudritzki et aT1ll989l) . The density and 
velocity distributions are related via the mass conservation equa- 
tion. 

The emergent spec trum is computed using CMF_FLUX 
dBusche & Hilliej|2005h . and results from a formal solution of 
the radiative transfer equation using the atmosphere model pro- 
vided by CMFGEN. A depth dependent microturbulent velocity, 
varying from 10 km s _1 in the photosphere to 200 km s _1 in the 
upper wind, is used. Proper pressure and Stark broadening is in- 
cluded for H, He i, Hen and CNO resonance lines. Other lines 
are treated by Gaussian profiles. 

For our purpose, we selected two sets of model: one with 
T eS = 30500 K and one with 42000 K. They corre spond to late 
and early O stars respectively dMartins e t al. 2005 : [). The atomic 
data for carbon were taken from various sources: the photoion- 
ization rates are from the opacity project dSeatonll987lPl oscilla- 
tor strengths and the low-temperature dielectronic transition data 
were kindly computed by Peter Storey (private communication). 
For collisional excitation we write the collisional excitation rate 
(in cgs units) as 



C(i,k) 



8.63 x 10- 8 ^^^ 



10000 



where N e is the electron density, g, the statistical weight of 
the lower level, the excitation energy, T the temperature, 
and Cl(i, k) the dimensionless collision strength. The collisional 
data for the six lowest l evels (singlet and triplet states) are from 
iBerrington et al. ( 1985b: for other colli sional data we use the ap- 
proximation of van Regemorterl d 19621) who writes the collision 
strength for a transition between two levels in-terms of the oscil- 
lator strength (for optically allowed transitions). 

3. C in 5696 

Cm 5696 (2s3p 1 P°-2s3d 'D) can either be in absorption or 
emission, and exhibits a complex variation with T e s and log g. 
To illustrate these dependencies, we show the C m 5696 profile 
for r eff = 30500 K and logg=3.25 and 3.75 in Fig.Q](left panel). 
In the low \ogg model Cm 5696 is in emission, while in the 
high log g model it is in absorption. The opposite behavior is ob- 
served in a much hotter model with = 42000 K - C m 5696 
is in absorption for logg = 3.5 but in emission for \ogg = 4.0. 

As C in 5696 is primarily photospheric a key quantity con- 
trolling the line strength (and whether the line is in emission or 
absorption) is the line source function. An increase in the line 
source function will weaken an absorption line, possibly driving 
the line into emission. Conversely, a decrease in the line source 
function will enhance (weaken) the absorption (emission) pro- 
file. The line source function can be increased by decreasing the 



1 Data for the opacity project can be obtained from 



http://cdsweb.u-strasbg.fr/topbase/testop/TheOP.html 
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Figure 1. Spectra of C m 5696 illustrating the dramatic dependence of the line profile on surface gravity and effective temperature. 
The left panel is for T e g =30500 K while the right panel is for 7^=42000 K. Particularly striking is the strong dependance of the 
profile on log g, and that this dependance is of opposite sign in models of different effective temperature. 



population of the lower level, and/or by increasing the popula- 
tion of the upper level. Thus, to understand the origin of C in 
5696, we need to understand the transitions populating and de- 
populating its lower and upper level. 
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Figure 2. Simplified Grotrian diagram for the singlet system of 
C m. The C ra 5696 transition is indicated in red. 



Fig. [2]presents a simplified Grotrian diagram of the singlet 
system of C m. We see that the lower and upper levels of C in 
5696 (2s3p 'P° and 2s3d 'D respectively) are connected to lower 
levels by several extreme UV transitions, and to upper levels by 



far UV/UV transitions. Level 2s3p 'P° is also connected to the 
ground level by a transition at 386 A. 

To illustrate which transitions control the lower and upper 
level populations of Cm 5696 we computed different models 
(all with r eff =30500K, logg=3.25, M = 2.0 x 10~ 6 M and 
Vco = 2000 km s ) in which we selectively removed specific 
C in transitions by reducing their oscillator strength by a factor 
10000. The results are displayed in Fig.[3](left panel). As an ex- 
ample, the dotted green line is the model in which C m 884 was 
removed. This figure highlights the key role of Cm 884. Its ef- 
fect is dramatic: C m 5696 switches from emission to absorption 
when it is removed from the model. Several other C m lines (e.g., 
Cm 386 and Cm 574) also have a significant impact on Cm 
5696, although with a reduced effect compared with C m 884. 

The right panel of Fig. [3] shows the departure coefficients 
of levels 2s3d 'D and 2s3p 'P (upper and lower levels of Cm 
5696) in the initial model and the models without C m 884 and 
Cm 574 (i.e., the lines leading to the strongest changes). We 
see that in the absence of the C in 884 A transition, the depar- 
ture coefficients are close to 1.0, while in the initial model, the 
lower level is depopulated compared to the upper level, leading 
to emission. The total radiative rate0 in C in 884 is positive in 
the initial model (see left panel of Fig. |5). Thus, C m 884 acts 
as a drain for level 2s3p P°. When we remove it, the drain is 
suppressed and the 2s3p P° population is higher, as seen in Fig. 
|3] (right panel). 

Similarly, the total radiative rate in C m 574 is also positive 
(see right panel of Fig. |5j, draining the population of 2s3d 'D. 
When C m 574 is removed, 2s3d 'D is less depopulated, its de- 
parture coefficient is higher than in the initial model (Fig.0 and 
C in 5696 emission is strengthened. 



2 The total radiative rate is the difference between the number of ra- 
diative transitions from the upper level towards the lower level and the 
number of radiative transitions from the lower level towards the upper 
level. It is positive when there is a net decay of electrons from the upper 
state and negative when there is a net excitation of electrons into the 
upper state. 
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Figure 3. Effect of the C m UV lines on the formation of C in 5696. Left: Initial C m 5696 spectrum (solid) and spectra from models 
in which several C in UV lines have been artificially removed. Right: Departure coefficient of the C m 5696 lower and upper levels 
in the initial model (solid lines), the model without C in 884 A (dashed lines) and the model without C m 574 A (dotted lines). 
Departure coefficients are plotted as a function of Rosseland opacity. The dot-dashed line at the bottom indicates the C in 5696 
formation region expressed in terms of the contribution function as a function of TR OSSe i an( j. 

lines on C m 5696 for the model with = 42000 K and log g = 
3.5. Compared to the left panel of Fig.El(Te ff = 30500 K, log g = 
3.25), lines have qualitatively similar effects. Removing C in 574 
leads to a stronger emission than the initial model; removing C m 
386 also leads to a stronger emission, but not as strong as in the 
case of Cm 574; and removing Cm 884 leads to more absorp- 
tion. But at 42000 K, it is the removal of C m 386 and C m 574 
(and not C m 884) which leads to the switch between absorption 
and emission. Thus depending on the temperature the C m 5696 
emission is driven by different C m UV lines. In appendix IB"1 we 
present the effect of these lines on the appearance of C m 8500 
(its upper level is the lower level of C m 5696) relevant for future 
Gaia observations. 



3.1. Explanation of log g effect 

Earlier we saw that the C m 5696 profile varies dramatically with 
both T e ff and log g. We can now use our understanding of the key 
role of C m 386, C m 574 and C m 884 in controlling the strength 
of C in 5696 to explain this behavior. 

In the left panel of Fig. [5] we show the total radiative rates 
from level 2s3p 1 P° (i.e. the lower lower level of Cm 5696) for 
Figure 4. Same as left panel of Fig.[3]for the model with r eff = logg = 3.25 and 3.75 at 30500 K. Focusing on the low gravity 
42000 K. case, we see that over the C in 5696 formation region the highest 

rate is observed for the C in 884 A transition. In addition, this 
rate is positive, indicating a depopulation of the 2s3p 'P level 
In the model with r eff =42000 K0 the physical conditions are as seen previously. Turning to the log g = 3.75 model, we see that 
different and the reaction of C m 5696 to the removal of UV lines the C m 884 A radiative rate has a qualitatively different behavior 
is slightly different. Fig. [4] shows the effect of different C m UV - in most of the C m 5696 formation region it is negative, and 

hence it does not act as a significant drain on the population of 

3 At r cff =42000 K we have adopted the same mass loss rate as in the the lower level of C m 5696. Consequently, C m 5696 switches 
low r c ff model for a better comparison. to absorption. 
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Figure 5. Effect of log g on the formation of Cm 5696 for r e ff = 30500 K. Left: Total radiative rates out of level 2s3p 'P at logg 
= 3.25 (upper panel) and logg = 3.75 (lower panel) as a function of Rosseland opacity. Notice how the largest rate is for A 884, 
and that this rate shows a qualitatively different behavior in the two models. The dotted lines indicate the C m 5696 line formation 
region. Right: Total radiative rates out of level 2s3d 'D at \ogg = 3.25 (upper panel) and \ogg = 3.75 (lower panel). 




We can understand the differences as follows. A decrease of 
log g translat es into a smaller Lyman ju mp and to a higher ioniz- 
ing flux (e.g. lAbbott & Hummeil l 19851). This also follows from 
density arguments — at a given T^, a low log g model is less 
dense, hence there are less recombinations and the atmosphere 
is more highly ionized than a higher logg model. Consequently, 
the ionization is stronger in the logg = 3.25 model atmosphere. 
C iv is the most abundant ion while at log g = 3.75, it is C m. In 



addition, 2s3p 1 P° and 2p 2 'D (the lower level of Cm 884) are 
collisionally coupled. 

The switch from absorption to emission for C m 5696 when 
\ogg decreases is due to the following chain of events: the 
ground state and 2p2p 1 D are depopulated due to the higher ion- 
ization and consequently the downward radiative rate in C m 884 
has increased; the Cm 884 drain on the 2s3p P° (lower level of 
C m 5696) population increases, which leads to emission in C m 
5696. Fig.|6]confirms that when log g decreases, the ground state 
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Figure 6. Effect of log g on departure coefficients. Departure co- 
efficients are plotted as a function of Rosseland opacity for the 
levels on C m 5696 (upper panel) and for the fundamental level 
(2s 2 l S) and the lower level of C m 884 (2p 2 *D). The dotted line 
shows the C in 5696 line formation region. 



and 2p2p D are depopulated, leading to a strong depopulation 
of 2s3p 'P through C m 884. 

At r eff = 42000 K, the effect of log gonCm 5696 is opposite 
(Fig-ID- The total radiative rates of the low and high log g models 
at 42000 K are shown in Fig. [7] C m 884 always acts as drain, 
and drives C m 5696 into emission. However, the resonance line 
C in 386 switches from near radiative equilibrium at \ogg = 4.0 
to a net negative rate at log g — 3.5 (left panel of Fig.|7]i, implying 
that it populates the 2s3p 'P° level. This compensates for the C m 
884 drain, and C in 5696 is in weak absorption. 

In conclusion, the reaction of C m 5696 to changes in log g 
is directly related to the line formation mechanism and the key 
role of the C m UV lines. As the temperature increases, the peak 
of the spectral energy distrubution shifts to shorter wavelengths 
and consequently higher energy C m UV lines become more and 
more important at the expense of lower energy transitions. 

3.2. Dielectronic recombinations 

So far we have focused on detailed non-LTE radiative transfer 
effects to explain the morphology of C in 5696. Dielectronic re- 
combination is another mechanism which might have an influ- 
ence on the appearance of C m 5696. As low temperature dielec- 
tronic recombination (LTDR) for levels permitted to autoionize 
in LS-coupling is handled via resonances in the photoionization 
cross-sections they cannot be easily omitted from the calcula- 
tions. Instead, we run a model in which we utilize the Einstein 
A values computed using intermediate coupling by Peter Storey 
(private communication) to add important LTDR transitions into 
the photoionization cross-sections. Although this means that we 
are exaggerating their influence (since we are double counting) 
we can still use the change in line profile to infer whether LTDR 
is important. Fig. [8] shows the effect of low temperature dielec- 



Figure 8. Effect of low temperature dielectronic recombination 
(LTDR) on Cm 5696 in the models with T eff = 30500 K (upper 
panel) and 42000 K (lower panel). The red solid curve shows the 
standard model profile, while the blue dashed profile shows the 
effect of "double counting" the influence of LTDR. The green 
dot-dashed line is the standard model in which the LTDR from 
levels forbidden to autoionize in LS coupling have been omitted. 



tronic recombinations on the appearance of C m 5696 at high and 
low effective temperature. The additional dielectronic recombi- 
nations, when included, lead to a slight increase of the emission 
line, but quantitatively, the effect is smaller than that of UV lines 
(see Figs.|3]and|4|i. We also run a test model in which we switch 
off LTDR for transitions not permitted to autoionize in LS cou- 
pling. These were also found to have only a small effect on C m 
5696. Thus LTDR is not a viable mechanism to explain the ob- 
served behavior of C m 5696 in O stars. Similar resul ts are found 
for N m 4634-4640 bv lRivero Gonzalez et all d2011l) . 

3.3. Line blanketing effects 

3.3.1. Microturbulence and metallic line interactions 

In Fig. [9] we show the effect of a change of the microturbulent ve- 
locity in the atmosphere model (and not in the calculation of the 
spectrum) on the C m 5696 line. A reduction of v tul b translates 
into a stronger emission. Modifying v tm b affects the efficiency 
of line-blocking and thus the temperature structure and EUV 
radiation field. But different microturbulent velocities also im- 
ply different metallic line width and thus different line overlap, 
which can affect radiative transfer in UV C m lines. 

We have seen in section [3] that the Cm 386, Cm 574 
and Cm 884 lines significantly influence the strength of Cm 
5696. Depending on the microturbulent velocity the lines have 
a typical width between 15 and 30 km s _1 , and hence will 
interact with other lines whose separations are less than 15- 
30km s . To test the effect of metallic lines on the Cm 
UV features and thus on Cm 5696, we have run test mod- 
els with T e ff = 30500 K. We selectively removed metallic fea- 
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Figure 9. Effect of the microturbulent velocity on the C m 5696 
emission strength in three CMFGEN models. The microturbu- 
lent velocity has been changed only in the atmosphere compu- 
tation. The spectra have all been computed with v tur b increasing 
from 10 to 200 km s throughout the wind. 



Figure 10. Effect of metallic lines on Cm 5696. The solid line 
is the initial spectrum. The dashed line is a model in which the 
oscillator strength of Fe v 386.262 A is artificially reduced by a 
factor 10000. The dot-dashed (dotted) line is the initial model in 
which Fe iv 574.232 A (S V 884.531 A) is reduced by the same 
factor. 



tures (Fe vi 386.186, Fe iv 386.21 1, Fe v 386.262, Fe iv 574.229, 
Fevi 574.230, Feiv 574.232, Ariv 574.243, Niiv 574.291, 
Nm 574.298, Nm 884.466, Feiv 884.497, Feiv 884.520, Sv 
884.530, and Fev 884.578) in the model atmosphere and com- 
puted the resulting C m 5696 line profile. Atomic data for the 
transitions found to influence C m 5696 and C in are discussed in 
AppendixlAl 

In all but three cases, which are shown in Fig. QJJ the 
C m 5696 feature is basically unchanged. The model with Fe iv 
574.232 removed (dot-dashed) has a stronger C in 5696 than the 
initial model (solid). The explanation for this emission increase 
can be inferred from Figs. QT] and [12] In Fig. QT| the departure 
coefficients of the upper and lower levels of C m 5696 are shown 
for both the initial model and the model without Feiv 574.232. 
Initially, the lower level is more depopulated than the upper 
level, explaining the emission. When Feiv 574.231 is removed, 
the upper level becomes more populated while the departure co- 
efficient of the lower level is barely changed. This translates into 
the stronger emission seen in Fig.fTUl 

The reason for the enhanced emission can be understood 
from Fig. Q~2] which shows the total radiative rates in C in 
574.281, Cm 5696 and Feiv 574.231. The total rate in Cm 
574.281 is positive in the initial model while the rate in the iron 
transition is negative. Consequently the Fe iv 574.232 transition 
removes photons from the C m 574.281 line. These photons are 
not available for re-absorption in the C m 574.28 1 and hence the 
population of 2s3d 1 D, the upper level of C m 5696, is reduced. 
When the Fe iv 574.232 line is artificially removed, it does not 
"steal" photons any more from the Cm 574.281 transition, re- 
sulting in a smaller depopulation of 2s3d 'D, as seen in Fig.fTTl 
and the emission in C m 5696 is stronger. 



The scenario presented above also partljQ explains why a 
lower microturbulent velocity results in a stronger Cm 5696 
emission. In the case of a large v tul -b, the line overlap between 
C m 574.28 1 and Fe iv 574.232 is wider, favoring the drain ef- 
fect described above. When v tm b is reduced, the reduction of the 
overlap translates in a less efficient drain of the 2s3d 'D level 
population, implying a stronger C in 5696 emission. 

In Fig. [10] we have also plotted the spectrum resulting from 
a model in which Fe v 386.262 A was reduced artificially. As 
for the Fe iv 574.232 line, the result is a stronger emission. This 
Fe v line can interact with C m 386 (Fig. [2). The total radiative 
rates in the initial model are negative for both lines (Fe v 386.262 
and C in 386). This means that the iron line absorbs photons that 
could in principle be absorbed by the carbon line. When the iron 
line is removed, the photons are effectively absorbed by the car- 
bon line. This excess of absorption could in principle overpop- 
ulate the 2s3p 'P° level and lead to a weaker emission. But in 
practice, the 2s3p P° population is reduced. This is explained 
by the need to maintain the ground level population constant, as 
indicated by our models (the ground level population being con- 
trolled by ionization/recombination processes). Since in the ab- 
sence of Fe v 386.262 the absorption in C m 386 is stronger and 
since the ground level population is unchanged, an increase of 
the downward transition is needed to maintain the ground level 
rate equation equilibrium. Consequently, the 2s3p P° is depop- 
ulated, leading to a stronger emission. Here again, this mech- 
anism explains that a smaller microturbulent velocity leads to 
less absorption by Fe v 386.262 and thus to a stronger C m 5696 
emission. 



4 The changes in the atmospheric structure - such as the temperature 
structure alteration - also play a role. 
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Figure 11. Effect of a change in the oscillator strength of Feiv 
574.232 A on the departure coefficient (plotted as a function of 
Rosseland opacity) of the levels of the C m 5696 transition. The 
red curves correspond to the upper level and the blue ones to the 
lower level. The solid lines are for the initial model. The dashed 
lines are the level populations in the model with a reduced Fe iv 
574.232 line. The dotted line is the C m 5696 formation region. 
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tion has been removed in the atmosphere calculation. The dotted 
line indicates the C in 5696 line formation region in the initial 
model. 



Finally, Fig. [TO] also shows the effect of the S V 884.531 
line on C m 5696. Removing this feature translates into a reduc- 
tion of the C in 5696 emission. The explanation comes from the 
Grotrian diagram in Fig. [2] and the effect of the C in 884 tran- 
sition on Cm 5696 (at T e ff = 30500 K). The total radiative rate 
in the S V line is negative throughout the C m 5696 line forma- 
tion region, while the C in 884 A line rate is positive. Thus the 
sulfur line enhances the drain of 2s3p 'P through C m 884. By 
removing the sulfur line, we reduce the drain, thus re-populate 
the 2s3p P° level, and consequently reduce the C m 5696 emis- 
sion. 

3.3.2. Metallicity effects 

We have seen in Section [37X1 that Cm 5696 was indirectly sen- 
sitive to the ionizing flux. We thus expect that any reduction of 
the Lyman continuum flux will affect the C m 5696 morphology. 
Since an increa se of the metal content sever ely reduces the ioniz- 
ing fluxes (e.g. ISchaerer & de Kotedfl997l) . and thus the ioniza- 
tion structure, we expect metallicity to have a significant impact 
on the appearance of C in 5696. 

Fig. [T3l shows the influence of the metallicity Z on the Cm 
5696 line profile and departure coefficients at r e ff = 30500 K. 
We have computed models with a global metallicity reduced by 
a factor five compared to the initial, solar metallicity model. We 
have also computed a model with a global metallicity of l/5' /! the 
solar metallicity, but with the carbon abundance fixed to the solar 
value. As expected, at low metallicity the ground state level is 
depopulated (see right panel of Fig.Qj). The collisional coupling 
with the 2p2p 'D level leads to a strong drain through C m 884 
A and to a depopulation of 2s3p 'P . But at the same time, we 
observe a strong depopulation of 2s3d 'D as well, so that the 
net outcome of a reduction of Z is a decrease of the C in 5696 
emission. This is not just an effect of the scaling of the emission 
strength with the carbon abundance, since in the model with a 
solar C/H but reduced metallicity we also observe a reduction 
of the line intensity. The reason is the depopulation of 2s3d 'D. 
It is due to a stronger radiative rate in the C m 574 line at low 
metallicity. This is explained by the reduced strength of the Fe iv 
574.232 line and consequently its reduced interaction with Cm 
574 as described in Sect. 13.3.11 

Appendix [C] shows how C m 5696 reacts to a change of 
metallicity for T e g = 42000 K. A combination of lower ionizing 
flux and reduced metallic line interaction explains the observed 
behavior. This highlights the complexity in predicting the shape 
of C in 5696. Metallicity can also affect the line formation by its 
influence on the wind strength (see Sect. P5J. 

3.3.3. Influence of N iv /1322 

lGauzid(fl966l ) argued that C m 5696 emission could be due to the 
interplay between Niv 322 and Cm 322 between levels 2s 2 ! S 
and 2p3s 'P°0, the latter being connected to 2s3d 'D (the upper 
level of C m 5696) by a transition at 2982 A (see discussion in 
iNussbaumerll 19711) . We have tested this possibility by removing 
C m 2982. The result, shown in Fig. [3] indicates that the 2982 A 
feature plays little role in populating 2s3d 'D. Hence, inducing 
the C in 322 transition and decay through C m 2982 by N iv 322 
is not responsible for the appearance of C m 5696 emission. 



5 For clarity and given the conclusion of our test, we have not shown 
these transitions in Fig. [2] 
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Figure 13. Effect of the metallicity Z on the formation of C in 5696 at T eB = 30500 K. Left: Spectra of C in 5696 at Z (solid), 1/5" 1 
Z (dashed) and l/5' /l Z but with C/H at the solar value (dotted). Right: Departure coefficient of the Cm 5696 lower and upper 
levels at Z Q (solid lines), l/5 f/l Z G (dashed lines) and 1/5''' Z with C/H at the solar value (dotted). The dot-dashed line indicates the 
line formation region. 



In conclusion, we have identified a general mechanism re- 
sponsible for the morphology of Cm 5696. It is the influence 
of UV C in transitions on the population of the upper and lower 
levels of C m 5696 that play a key role. The main lines are C m 
386, Cm 574, and Cm 884. At low r eff , Cm 884 is crucial for 
determining the lower level population, while at high , C m 
386 is more important. Subtle line-blanketing effects modify the 
strength of the C m 5696 emission. Obtaining an excellent fit to 
C m 5696 is an extremely challenging task since it requires the 
knowledge of several stellar parameters including the microtur- 
bulent velocity, the metallicity, as well as accurate wavelengths 
and oscillator strength of several weak metal lines. 



4. C in 4647-50-51 

We now turn to the study of Cm 4650. We follow the same 
method as that used to identify the formation mechanism of C m 
5696. 



4.1. Line formation mechanism 

The C in 4650 lines are produced by the triplet system of C m. 
They are the triplet equivalent of the C m 8500 transition in the 
singlet system (i.e., a transition between levels 3s and 3p). The 
triplet equivalent of Cm 5696 is Cm 9710. Fig. [141 shows the 
simplified Grotrian diagram for this system. 



3x10 s 



2xl0 5 



10 5 




Figure 14. Simplified Grotrian diagram for the triplet system of 
Cm. The Cm 4647-50-51 transitions are shown in red. The star 
symbols indicate multiplets. 



4.1.1. Effects of UV lines 

To understand the formation of C m 4650, we proceeded as for 
C m 5696, looking for key transitions that control the upper and 
lower level populations. Fig. [TBI shows the Cm 4647-50-51 pro- 
files in the low (upper panel) and high (lower panel) T e jf case. 
The solid line is the initial model, while the other lines corre- 



spond to models in which specific C m UV lines connected to 
either the upper of the lower level of C m 4650 have been re- 
moved (by reducing their oscillator strength by a factor 10000). 
Only the models with the largest changes are shown. The fol- 
lowing lines were tested: C m 538, Cm 818, C m 981, C m 1 175, 
Cm 1256, Cm 1428, Cm 1578, Cm 1620, Cm 1908, C m 2009. 
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Figure 15. Effect of Cm UV lines on Cm 4647-50-51 in the 
model with T eff = 30500 K (upper panel) and 42000 K (lower 
panel). The solid line is the initial model. The other lines are 
models in which a given C m UV line has been removed by re- 
ducing its oscillator strength by a factor 10000. 

The C in 538 A transition is found to play a crucial in the 
formation of the Cm 4647-50-51 line, since removing it leads 
to a strong increase of the absorption. The C m 2009 line also 
has a non negligible role on the shape of C m 4650, but its im- 
pact is much weaker than that of C in 538. The same results are 
obtained in the model with r e ff = 42000 K - all models have 
Cm 4650 in emission, except the model in which Cm 538 is 
removed and where absorption is observed. From the Grotrian 
diagram (Fig. RBI) we can suspect C m 538 to be a drain for level 
2s3s 3 S. Inspection of the departure coefficient of 2s3s 3 S in the 
initial model, and the model without C in 538 (not shown), con- 
firm that 2s3s 3 S is more populated when C m 538 is absent. 

4.1.2. Effect of log g 

This mechanism (depopulation of 2s3s 3 S by C m 538) explains 
the behavior seen in Fig. [16] The effect of log g on C m 4647-50- 
5 1 in a model at 42000 K is shown. The reduction of the gravity 
turns the weak emission in the log g = 4.0 model into a strong 
emission in the \ogg = 3.5 model. The departure coefficients 
are shown in the right panel of Fig. [16] While the upper level 
has b, close to 1 in both models, the lower level has a departure 
coefficient much lower than 1 in the low gravity model. This 
increases the line source function, leading to strong emission in 
the low gravity model. 

Given the Grotrian diagram seen in Fig.[14]one can conclude 
that the mechanism is similar to that of C m 884 for C in 5696. 
The higher ionization of the low log g model leads to a depopu- 
lation of the ground level (2s 2 'S) and the 2s2p 3 P° to which it is 
collisional coupled. The lower population of 2s2p 3 P° enhances 
the C in 538 drain (the total radiative rate is positive) depopulat- 
ing 2s3s 3 S and driving Cm 4647-50-51 into emission. Hence, 
C m 538 controls the population of the lower level of C m 4650 



and its morphology (absorption or emission) at both high and 
low effective temperature. The situation for C m 4650 is clearer 
than for C m 5696, since the only line affecting significantly its 
morphology is Cm 538. In particular, Cm 818, connected to 
2s3p 3 P° (upper level of C m 4650), has no impact on the ap- 
pearance of C m 4650 (Fig. IT~5T >. The reason is related to a quirk 
of the atomic physics: the Einstein coefficient A for the 818 tran- 
sition is only 2.3 x 10 6 compared with 7.2 x 10 7 for Cm 4650 
and 4.1 x 10 8 for Cm 538. 

As for C in 5696, we also tested the effect of low tempera- 
ture dielectronic recombination and found out that its role was 
marginal. The line profiles show slightly more emission (or less 
absorption) when LTDR is included, but the changes are mi- 
nor compared to the effects of UV lines highlighted above. In 
O stars, LTDR is not responsible for the strength of C m 4650. 

4.2. Line blanketing effects 

4.2.1. Microturbulence and metallic lines 

As for C m 5696 we find that microturbulence and line blanket- 
ing effects influence the morphology of C m 4650. In particular, 
an increase of v tur b lead to a significant strengthening of the ab- 
sorption. Inspection of our line list also indicates that several 
metallic lines are located near the C m 538 transition - the only 
UV line that significantly influences the strength of C m 4650. 

We tested the influence of the following lines on the strength 
of Cm 4650: Feiv 538.044, Feiv 538.057, Fevi 538.097, 
Fe vi 538.1 12, Fe iv 538. 176, Fe iv 538.265, Fe vi 538.285, Om 
538.320, Arm 538.348. Only Feiv 538.057 leads to a change 
of the C in 4650 profile - its absence translates into a "small" 
increase in the absorption (Fig. [TTT l. While the effect of Fe iv 
538.057 is important Cm 538 has a much larger impact. The 
Fe iv 538.057 has a negative total radiative rate and consequently 
assists in the depopulation of the 2s3s 3 S state, like Fe iv 574.232 
for the upper level of C m 5696. 

4.2.2. Metallicity effects 

Fig. Q~8] shows the effect of metallicity on C m 4650 at T e s = 
30500 K (left) and 42000 K (right). At low effective temperature, 
a global reduction of the metallicity translates into a stronger 
absorption. If the carbon content is kept at the solar value, the 
absorption is further strengthened, simply because there is more 
carbon and thus more C m in the atmosphere. The increase of the 
absorption due the reduction of Z can be understood as follows. 
We have seen in Fig.[17]that removing C m 538 and Fe iv 538.057 
lead to an increase of the C m 4650 absorption. A reduction of Z 
corresponds to a reduction of the strength of C m 538 and Fe iv 
538.057. Consequently, a lower Z implies a stronger Cm 4650 
absorption. 

For T eff = 42000 K, the behavior of C m 4650 is the same 
as that of Cm 5696 (see Appendix [Cb. Reducing Z produces 
a weaker emission. But keeping C/H at the solar value and re- 
ducing the abundance of all other elements lead to a C m 4650 
emission stronger than at solar metallicity. The reason is that the 
condition for emission are more favorable at low Z. In particu- 
lar, the ionization is higher, draining the population of 2s3s 3 S 
significantly even if C m 538 is weaker at low metallicity. 

As highlighted in the introduction, C m 4650 and C m 5696 
are usually the strongest C m lines in O stars optical spectra. It is 
tempting to use them to perform carbon abundance determina- 
tions. However, we have seen that both lines critically depends 
on metallic lines the characteristics of which are not accurately 
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Figure 16. Effect of log g on C m 4647-50-5 1 . Left panel: the red line is a model with log g = 3.5; the blue line is the same model 
with log g = 4.0. The effective temperature is 42000 K. Right panel: departure coefficient of the upper and lower levels of C m 
4647-50-51 (upper panel) and of levels 2s 2 'S (n=l), 2s2p 3 P° (n=2,3,4), 2p2p 3 D (n=6,7,8) (lower panel) for logg = 3.5 (solid 
line) and \ogg = 4.0 (dashed line), plotted as a function of Rosseland opacity. The dotted line is the C in 4650 formation region. 




Figure 17. Synthetic spectra around the C m 4650 triplet in mod- 
els with T e g = 30500K. The red line is the initial model. The blue 
line is the same model in which we have removed the C m 538 A 
line. The green model is the same as the initial model in which 
the Fe iv 538.057 A line has been removed. 



Figure 18. Effect of metallicity on the C m 4650 triplet at 30500 
K (left) and 42000 K (right). The solid line is a solar metallicity 
model. The dashed line is a model with Z = 1 /5 r/, Z . The dotted 
line is the same model as the dashed line model but with a solar 
C/H. 



known. We illustrate in Fig.[19]how such uncertainties can affect 
carbon abundance determinations. Removing Feiv 574.232 in a 
model with C/H = 1/2 C/H Q leads to a C m 5696 profile almost 
as strong as in the solar C/H model (left panel). For Cm 4650, 



the effect of removing Fe iv 538.057 is even larger since the ab- 
sorption in the C/H =1/2 C/H Q without this iron line is stronger 
than the absorption in the solar C/H model. In conclusion, us- 
ing C in 4650 and C in 5696 to determine the carbon content of 
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the emission is seen. This is what is commonly observed in wind 
sensitive lines such as Ho-. In the 42000 K model, below 10~ 7 
M© yr -1 , Cm 4650 has a rather constant emission level, with 
slight fluctuations. Above 10~ 7 M© yr _1 , the emission increases 
with mass loss rate. 

A detailed explanation of the observed behavior is provided 
in Appendix [D] Basically, it is related to the influence of UV 
lines on the populations of the upper levels. As we increase 
the mass-loss rate we increase the velocity at a fixed Rosseland 
optical depth, and this alters the radiation transfer in the UV 
lines. For example, at higher velocities the lines desaturate ear- 
lier which can enhance both photon escape and photon pumping. 
Which of the two processes dominates depends on both the radi- 
ation field in the lines neighborhood and the line under consider- 
ation. These two processes have opposite effects on the strength 
of an optical line - photon escape in the UV line enhances optical 
emission (weakens absorption) while photon pumping enhances 
absorption (weakens emission) when the UV line's upper state 
is the lower level of the optical transition. The converse occurs 
when the UV lines is connected to the upper state of the optical 
transition. The behavior of C ni 5696 at low mass loss rates is 
complex and not straightforward. 



Figure 19. Effect of C/H on the intensity of C m 4650 (left) and 
C m 5696 (right) at T eff = 30500 K. In both panels, the solid line 
is a solar metallicity model. The dashed (dotted) line is a model 
with C/H = 1/2(1/5) C/H , all other elements having the solar 
abundances. In the left (right) panel, the dot-dashed line is the 
C/H = 1/2 C/H model in which the Fe iv 538.057 (574.232) line 
has been removed. 



O stars is very risky, and quoting an uncertainty of less than a 
factor of two in the derived C/H using these lines is unrealistic. 

5. Wind effects 

We devote a special section to the effects of winds on the 
morphology of the Cm optical lines. Details are also given in 
Appendix iDl In Fig. [20] we show the effects of changes in the 
mass loss-rate at 30500 K and 42000 K. The effects are complex, 
sometimes counter-intuitive, and do not change monotonically 
with mass-loss rate. The changes show that the optical lines can- 
not be modeled without modeling the wind, and that an accurate 
estimate of the mass-loss is required. Because of clumping, ob- 
taining an accurate estimate of the mass-loss rate is not easy. 
Because mass-loss rates strongly correlate with luminosity, we 
expect the C m optical lines (at fixed T e ff and log g) to correlate 
with luminosity. 

Let us first examine C in 5696. At T e g = 30500 K, starting 
from the low mass loss rate model (3 x 10~ 9 M© yr _1 ), we 
see that the emission line profile first strengthens when M in- 
creases. A further increase of M from 10~ 8 to 3 x 10~ 8 M© 
yr~' translates into a change from emission to weak absorp- 
tion. Above 3 x 10~ 8 M© yr , the increase of M resumes to 
produce stronger emission. Thus, the behavior of the C m 5696 
line strength/morphology with increasing mass loss rate is not 
monotonic. At higher effective temperature, there is a clear anti- 
correlation of strength of C m 5696 emission with mass loss rate 
above 10~ 7 M© yr . Below that value, the emission stregnth is 
almost at the level of the 10~ 7 M© yr~' model. 

For C m 4650 at low T e ff , an increase of the mass loss rate is 
associated with a decrease of the absorption and/or an increase of 



6. Concluding remarks 

We have investigated the formation processes of C m 5696 and 
Cm 4650 in the atmosphere of O stars. We have used the non- 
LTE atmosphere code CMFGEN to perform simulations at two 
effective temperatures (30500 K and 42000 K) corresponding to 
late and early spectral types. Our results can be summarized as 
follows: 

• For wind densities typical of O stars, the formation of C m 
5696 is controlled by several UV transitions: Cm 386, Cm 
574 and Cm 884. The latter is crucial at low T e s: when it is 
removed from models, C m 5696 switches from emission to 
absorption. Cm 574 and Cm 386 play the same role at high 

• Cm 884 acts as a drain on the 2s3p 'P° population (lower 
level of Cm 5696). When it is removed, 2s3p l P° is 
more populated and emission is reduced/suppressed. Cm 
574 drains the 2s3d 'D population and Cm 386 populates 
2s3p 'P°. When they are removed, 2s3d 'D is more popu- 
lated and 2s3p 'P is depopulated, favoring emission. 

• Lower gravity models have larger ionizing flux. As a conse- 
quence, the lower level of C m 884 (2p2p 'F ), which is col- 
lisionally coupled to the ground level, is depopulated at low 
logg. The drain of 2s3p 'P by Cm 884 is thus increased. 
This leads to an increase of C m 5696 emission when log g 
decreases at 30500 K. At T eff = 42000 K, the effect of log g 
is opposite and is due to the role of Cm 386 in populating 
the lower level of C m 5696. 

• The strength of C m 5696 is very sensitive to the action of a 
few metallic lines on the radiative transfer near 386, 574 and 
884 A. In particular, Fe iv 386.262 and Fe iv 574.232, when 
removed from models, lead to a stronger C m 5696 emission. 
If S v 884.53 1 is removed, the C m 5696 emission is reduced. 
The reason is the overlap of these features with the key C m 
lines controlling the shape of C m 5696. As an example, Fe iv 
574.232 extracts photons from C m 574 leading to a lower 
population of 2s3d 'D. Consequently, Feiv 574.232 limits 
the C in 5696 emission. When it is removed, 2s3d 'D is more 
populated and the C m 5696 emission is stronger. 
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Figure 20. Effects of a change in the mass loss rate on C m 4650 and C m 5696 in models with T e g = 30500 K (left) and T e g =42000 
K (right). The mass loss rates (in Mq yr _1 ) are indicated in the C m 5696 panels. 



• The formation of C m 4650 is controlled by radiative transfer 
at 538 A. The C m 538 line drains the population of 2s3s 3 S, 
the lower level of Cm 4650. When Cm 538 is removed, 
2s3s 3 S is more populated and C in 4650 switches from emis- 
sion to absorption, or shows a stronger absorption. 

• Feiv 538.057 enhances photon escape in the Cm 538 line 
which helps to reduce the population of 2s3s 3 S, and thus 
weakens the Cm 4650 absorption (or enhances the emis- 
sion). 

• Low temperature dielectronic recombinations have a limited 
effect on the appearance of C m 5696 and C m 4650 in O 
stars. 

• When the mass loss rate increases, a growing fraction of the 
C in 4650 and C m 5696 lines is formed in the wind, leading 
to a larger emission. In addition, as M gets higher, the forma- 
tion region of the C m lines (UV and optical) shifts to higher 
velocities. The associated Doppler shifts lead to desaturation 
of the key C m UV lines. This affects the level populations 
and consequently the intrinsic C in 4650 and C in 5696 line 
profiles. At low mass-loss rates, the change in line profile 
with a change in the mass-loss rate is not monotonic. 

• Metallicity effects are a complex interplay between C m UV 
lines, metallic Fe iv and S v lines, ionization, winds and 
abundance effects. There is no simple way of predicting how 
C in 4650 and C m 5696 will react to a change of metallicity 
or abundance of specific elements. 

In view of our results, it seems highly risky to derive accurate 
carbon abundances from C in 4650 and C in 5696 in O star atmo- 
spheres. Uncertainties of at least a factor of two on C/H are ex- 
pected just because of the uncertainty on atomic data and wave- 
lengths of metallic lines. Priority should be given to other carbon 
lines, such as the Cm 2s4f 3 F°-2s5g 3 G multiplet at 4070 A, or 
even UV lines (Civ 1164 and Cm 1176) provided the wind is 
not too strong and they are mainly photospheric. 

The current results can be used to try to understand the be- 
havior of C m 4650 in Of?p stars, many (all?) of which are mag- 



netic (iDonati et alj2006trHubrig et alj2008tlMartins et al.l2010t 
IWade et aUl201 1 al ibi). Of?p stars show periodic variability, espe- 
cially in Cm 4650 emission which can be as strong as in the 
neighbouring N m lines before almost vanishing. H ff shows the 
same variability a nd is strongest when C m is at its maximum 
emission strength. ISundqvist et al.l d2012l) demonstrated that the 
changes in H ff emission of HD 191612 could be understood as 
a result of viewing the dynamic magnetosphere from different 
angles. 

Due to wind confinement by the magnetic field, an overdense 
region is created in the magnetic equator. If the star is seen from 
the magnetic pole, a large surface of the overdense region is ob- 
served and the H ff emission is stronger than if the star is seen 
from the magnetic equator. In the present study, we have seen 
that C in 4650 can show emission if the wind contribution be- 
comes dominant, and/or if for some reason the drain effect of 
Cm 538 is increased so that 2s3s 3 S is depopulated. If we ob- 
serve HD 191612 from the magnetic pole, the wind velocity 
at a given radius is larger than in the magnetic equator, due to 
confinement (see Fig. 2 of Sundqvist et al. 2012). The role of 
the wind is thus stronger, favoring C in 4650 emission. In addi- 
tion, since the atmosphere density is lower in the magnetic pole 
direction, we can expect a larger ionizing flux due to reduced 
opacity. As seen in our study, this favors the drain of 2s3s 3 S 
through C m 538, and thus C m 4650 emission. Hence, the exis- 
tence of an equatorial overdensity could lead to the observed C m 
4650: more emission when the star is seen pole-on, less emis- 
sion when it is seen equator-on, just as H a . Of course, this needs 
to be tested by detailed 2D radiative transfer calculation, but the 
present results provide a good working frame. 
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Table A.l. Elements and number of levels (super-levels) in- 
cluded in the CMFGEN models. 





T cff = 30500 K 


r cff = 42000 K 


Hi 


30 (30) 


30 (30) 


He i 


69 (69) 


69 (69) 


Hi 


30 (30) 


30 (30) 


Cn 


39 (21) 


- 


Cm 


243 (99) 


243 (99) 


Civ 


64 (64) 


64 (64) 


Nn 


105 (59) 


- 


Nm 


287 (57) 


287 (57) 


Niv 


70 (44) 


70 (44) 


Nv 


49 (41) 


49 (41) 


On 


274 (155) 


- 


Om 


104 (36) 


104 (36) 


Oiv 


64 (30) 


64 (30) 


Ov 


56 (32) 


56 (32) 


Nen 


48 (14) 


48 (14) 


Nem 


71 (23) 


71 (23) 


Neiv 


52 (17) 


52(17) 


Sim 


50 (50) 


50 (50) 


Si iv 


66 (66) 


66 (66) 


S m 


78 (39) 


78 (39) 


S iv 


108 (40) 


108 (40) 


Sv 


144 (37) 


144 (37) 


Arm 


138 (24) 


346 (32) 


Anv 


102 (30) 


382 (50) 


Arv 


29 (14) 


376 (64) 


Cam 


88 (28) 


232 (44) 


Caiv 


72(19) 


378 (43) 


Fern 


607 (65) 


607 (65) 


Feiv 


1000 (100) 


1000 (100) 


Fe v 


1000 (139) 


1000 (139) 


Fevi 


1000 (59) 


1000 (59) 


Fe vn 




245 (40) 


Nim 


150 (24) 


150 (24) 


Niiv 


200 (36) 


200 (36) 


Niv 


183 (46) 


183 (46) 


Ni vi 


182 (40) 


182 (40) 


Ni vii 




67 (20) 



Appendix A: Metal line atomic data 

Table lA~n provides the number of levels and super-levels used 
in our computations. Basic atomic data for the important tran- 
sitions that influence C m 5696 and 4650 are listed in Table [A] 
There is some uncertainty in the oscillator strengths, and this 
must be regarded as a source of uncertainty when modeling the 
Cm lines. For a comparison we provide an alternative f values 
that is available in the literature for each transition. Fortunately, 
for the transitions which significantly influence C m 5696 and 
4650, the wavelengths and/or energy levels for the relevant tran- 
sition are known. The wavelengths of most of the other transi- 
tions tested are (very) uncertain, but potentially could have influ- 
enced the modeling because of line overlap. A more worrisome 
concern is there may be some other important transition which 
overlaps with one of the key UV transitions, and which has a 
wrong wavelength (because it has not been measured in the lab- 
oratory) or was omitted from the model atoms used in the model 
atmosphere. The wo rk presented in the current paper, and that of 
iNaiarro et al] d2006l) . highlights the need for accurate Fe wave- 
lengths and oscillator strengths for modeling O star spectra. 



Appendix B: Effect of metallic lines on C in 8500 

The Grotrian diagram for the C in singlets (Fig. [2]i shows that 
the C in 8500 transition should also strongly depend on the ef- 
ficiency of the Cm 884 and Cm 386 drains since 2s3p 1 P° is 
the upper level of all three transitions. Consequently, any mod- 
ification of the 2s3p L P° population will change the Cm 8500 
morphology. 

In Fig. IB.ll we show how this line reacts to the C m UV lines 
at r e ff = 30500 K (upper panel) and r eff = 42000 K (lower panel). 
As expected, removing the drain from C m 884 translates into a 
weaker absorption at low T e g and to emission at high T e g . C m 
386 populates 2s3p 'P at high 7^ and thus removing it leads 
to absorption in C in 8500. C m 690 has also a strong drain ef- 
fect on 2s3s 'S since removing it leads to a strong absorption. 
The C in 8500 transition is of special interest in the context of 
future observations of O stars by Gaia, since it is basically the 
only non H/He line in the spectral window of the Radial Velocity 
Spectrometer onboard the satellite. Understanding its properties 
is thus important to classify O stars observed by Gaia. 
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Table A.2. Basic line data for the important overlapping transitions 



Species 


Transition 


Wavelength(A) 


f (Model) 1 


f(Alt) 2 


Feiv 


3d 5 4 F 5/2 -3d 4 ( 3 D)4p 4 D° /2 


574.232 


0.00523 


0.0131 


Fe v 


3d 4 5 D„-3d 3 ( 4 F)4p 5 D° 


386.262 


0.0406 


0.0341 


Sv 


3s3d 3 D 3 -3p3d 3 D° 


884.531 


0.0266 


0.00745 


Fe vi 


3d 5 4 G 7/2 - 3d 4 ( 3 H)4p 4 G° /2 


538.057 


0.121 


0.0174 



1. Model oscillator strengths are from R.L. Kurucz and B. Bell (1995 Atomic Line Data) Kurucz C D-ROM No. 23. Cambridge, Mass.: Smithsonian 
Astrophysical Observatory. Newer data are available at http://kurucz.harvar d.edu/ dKuruczll2009h 

2. Alternate f values for Fe are from lBecker & Butler! ( 1995b. 1992, 1995a) while that for S v is from[Ralchenko et al. (201 2 (Version!) . 
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Figure B.l. Effect of Cm UV lines on the appearance of Cm 
8500 for r eff = 30500 K (top panel) and 42000 K (lower panel). 



Figure C.l. Effect of metallicity on C m 5696 in the model with 
r eff = 42000 K. 



Appendix C: Effect of metallicity on C in 5696 for T eS 
= 42000 K 

Fig. IC. 1 1 presents the effect of a change of metallicity on Cm 
5696 in the model with T e g = 42000 K. We show a model with a 
global reduction of Z by a factor 5 (dashed line) and a model with 
the same reduction, but in which C/H is kept at the solar value 
(dotted line). The solid line is the solar metallicity model. A 
global reduction of Z (including a reduction of C/H) by a factor 5 
leads to a significant reduction of the emission in C m 5696. But 
a model with Z = 1/5 rA Z and a solar carbon content has a C in 
5696 line slightly stronger than the solar metallicity model. The 
higher ionizing flux at low Z implies a depopulation of 2s3p P° 
and thus more emission. In addition, the reduced effect of Fe iv 
574.232, Fe v 386.262 and S v 884.531 tend to increase emission 
(see Fig. [Toll. All in all, the conditions are more favorable at low 
Z, explaining that for a similar C/H, the emission is stronger than 
at solar metallicity. 



Appendix D: An explanation for the influence of the 
wind on optical C in 5696 lines 

In this section we provide an investigation of the mechanisms 
responsible for the trend observed in Fig. [20] which shows the 
effect of mass loss on C in 4650 and C m 5696. 

We first focus on Cm 5696. Its beh avior in the 42000 K 
model can be understood as follows. Fig. ID. II shows the depar- 
ture coefficients of the lower and upper levels of C m 5696 in a 
model with M = 4x 10~ 7 M G yr _1 (dashed line) and 2x 10~ 6 M© 
yr~' (solid line). In the low M model, 88% of the line is formed 
at log TRosseiand > —1.8, corresponding to velocities smaller than 
10 km s . In that region, the lower level is always less popu- 
lated than the upper level, a favorable condition for emission. In 
the high M model, a significant fraction of Cm 5696 (41%) is 
formed at velocities higher than 10km s _1 (logTR OSSe i an d < -1.2). 
In the low velocity part of the formation region, the lower level 
is less populated than the upper level, but above 10 km s _1 , the 
lower level is much more populated. As a result, the global C m 
5696 emission is weaker than in the low M model. The behavior 
of C in 5696 with mass loss observed in the left panel of Fig. [20] 
is thus a direct consequence of the strengthening of the wind. 

The behavior of 2s3p 'P° and 2s3d 'D when the wind sets in 
is explained as follows. Fig. ID. 21 shows the total radiative rates 
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Figure D.l. Effects of mass loss rate on 2s3p 'P (blue) and 
2s3d 'D (red) departure coefficients in a model with = 42000 
K. The dashed (solid) lines corresponds to a model with M = 
4 x 10~ 7 (2 x 10~ 6 ) M yr _1 . The dot-dashed (short dash - 
long dashed) line is the line formation region in the low (high) 
M model. 



results from this, with the consequence of stronger C m 5696 
emission. All in all, a higher mass loss rate leads to a stronger 
C in 5696 emission at low T s g . 

Wind effects can also explain the behavior of the C in 4650 
triplet in Fig. [20] The higher the mass loss rate, the larger the 
fraction of the line formed at velocities higher than 10 km s~'. 
The wind contribution to the line profile increases, implying 
more emission at high effective temperature and less absorp- 
tion at low temperature. In addition, at high r e ff, the UV lines 
populating the upper and lower level of C in 4650 are formed at 
velocities larger in the high M model than in the low M model. 
They all have positive radiative rates over the line formation re- 
gion. As a result C m 538 depopulates 2s3s 3 S as soon as Doppler 
shifts allow the line to absorb more continuum photons. In the 
same time, Cm 1256, Cm 1428, Cm 1578, Cm 2009 populate 
2s3p 3 P°. The effect is an increase of the Cm 4650 emission, as 
observed in Fig. [20] 



in lines connecting the upper and lower levels of C m 5696. The 
rates are shown for two mass loss rates discussed above (4x 10~ 7 
and 2 x 10~ 6 Mq yr~'). When M is increased, we see 1) that the 
rate in the C m 574 line increases slightly (compared to the rate 
in C in 5696) and 2) the rate in C m 386 shows a slightly more 
negative excursion at the end of the C m 5696 line formation re- 
gion. The net effect is to re-populate 2s3p 'P° and to depopulate 
2s3d 'D, leading to a reduced emission. The higher rates in C m 
386 and C m 574 can be understood as a desaturation effect. Due 
to the higher mass loss rate, the formation region is shifted to- 
wards larger velocities. The associated Doppler shift allows Cm 
386 and C m 574 to absorb more photons from the adjacent con- 
tinuum, and consequently the radiative rates in those lines are 
increased. As we have seen in section 14.11 C m 884 does not 
play an important role in the formation of C m 4650 in the high 
temperature model. 

In the 30500 K model, the increase of the emission when M 
strengthens from 3 x 10~ 8 to 2 x 10~ 6 Mq yr~' is mainly due to a 
growing contribution of the wind to the line profile, as at higher 
temperature. The fraction of the line formed above 10 km s~' 
grows from 0.6% for M = 3 x 10~ 8 M© yr _ 1 to 23.5% for M = 
2 x 10~ 6 Mq yr~'. Consequently, the contribution of the wind 
(emission^) is larger for higher mass loss rates. In addition to this 
dominant effect, there is also a higher radiative rate in the C m 
884 line in the photosphere (below 10 km s~'). In the same time, 
the rates in the C m 386 and C m 574 lines are not significantly 
modified compared to the other rates. Desaturation is responsible 
for the enhanced rate in Cm 884. A lower 2s3p 'P population 

6 The wind contribution is the radiation coming from material above 
the photosphere. Radiation is emitted in all directions in spectral lines 
and a fraction of the photons are parallel to the sight line, thus adding 
to the light directly coming from the photosphere. 
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